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A  Statistical  Model  for  Hydroaen  Halide  Pro; 

Using  Inforiuation  Theory* 


"lout  ions 


D.  H,  Stone  and  R.  L.  Kerber 
Michigan  State  University 
East  Lansing,  MI  48824 


-Chemical  laser  modeling  is  dependent  on  rhe  reaction  rate 
coefficients  available  from  both  experiment  and  theory.  A 
statistical  model  has  been  developed  to  correlate  the  relative 
rate  coefficients  for  the  laser  pumpinc  reactiens: 

(I)  F  +  -  HF{v,J)  ^  H,  (II)  F  H-  D2  -  DF(v,J)  +  D, 

(III)  H  -r  F^  ^  HF(V,J}  r  F,  (IV)  D  +  F^  -  DF;v,J)  F, 

(V)  H  +  CL-  -  HCi.(v,J)  t  CL,  (VI)  D  +  Cl-  -  DCL(v,J)  +  cU 
(VII)  H  +  Br'^  HBriv,J)  ->-  Br,  and  (VIII)  D  -  Dr'^'  ^  DBr(v,J)  t  5r. 
The  detailed  product  distributions  for  Reactions  ^IV)  and  (VIII) 
are  generated  by  the  model*  these  distributions  ha'/e  not  yet 
been  experimentally  determined.  'The  model  uses  surprisal 
analysis  to  relate  the  product  rotational  distributions  for 
each  reaction  by  considering  each  v'ibrational  level  separately. 
Using  Polanyi's  exterimental  data  with  RRiiO  prior  rates  results 
in  rotational  surprisals  which  are  approximately  quadratic  in 
form  and  vary  in  width  with  vibrational  level.  .A  .model  was 
developed  and  applied  to  Reactions  (III)  and  (VII)  which  assumes 
a  reaction  complex  interaction  among  the  product  vibrational 
levels.  An  adjusted  number  of  product  states  made  available 
to  each  level  by  the  interaction  is  then  computed.  The  logarithm 
of  this  number  is  related  to  the  width  (or  entropy)  of  each 
rotational  surprisal  distribution  and  a  correlation  is  observed. 
The  model  also  predicts  the  degree  to  which  ihe  surprisals  s)<ew 
toward  high  rotational  levels,  ^he  model  results  coupled  with 
the  observed  vibrational  distributions  favorably  reproduce  the 
rotational  distributions  for  Reactions  (III)  and  (VII).  Assuming 
an  isotopic  independence  for  some  parameters  between  Reactions 
(III)  and  (IV)  and  between  (VII)  and  (VIII) ,  the  model  can 
generate  the  full  vibrotational  distributions  for  (IV)  ,ind  (VIII) 
from  a  small  set  of  input  parameters. 

*This  worl;  was  supported  by  AFOSR  Grant  No.  7  5-284  2  and  NSF 
Grant  No.  ENG  76~Oh'733. 


I.  INTRODUCTION 


Chemical  laser  modeling  is  dependent  on  the  reaction  rate 
coefficients  available  from  both  experiment  and  theory.  A 
comprehensive  computer  model  must  incorporate  potentially 
hundreds  of  rates  for  the  various  pumping  and  relaxation 
mechanisms,  in  order  to  accurately  predict  laser  performance. 
Available  reaction  rate  data  for  HF  chemical  lasers  is  taken 
from  selected  experiments  and  trajectory  calculations  as  re¬ 
viewed  by  Cohen  and  Sort  in  references  i  and  2.  Not  all 
reaction  rates  of  interest  in  the  HF  laser  have  been  studied 
and  signif  it-c»nU  ancei  r a j.n L-ies  are  present  in  many  that  are  known. 
Techniques  are  needed  to  expand  the  data  base  from  a  few 


accurately 

measured  re 

accion  rates  to  a 

complete 

rate  set. 

In 

this  paper. 

we  apply  t 

he  information  - 

theoretic 

or  "surpris 

al" 

approach  to 

reaction  p 

rcauct  distribut i 

ons  as  dev 

■eloped  by 

Bernstein , 

Levine,  and 

ben-bhaul  (references  3-5) 

to  the 

exper imenta 

1  distribut 

ions  obtained  by 

Polanyi,  b 

.’oodall,  and 

Sloan  (refe 

rences  6-7) 

for  the  pumping 

reactions 

F  +  H., 

-  HF(v,J)  +  K 

(I) 

F  t 

-  DF(v,J)  +  D 

(ID 

H  +  F^ 

*  HF(v,J)  +  F 

(III) 

Our  objectives  are  to  develop  a  model  which  correlates  tne 
vibrotational  product  .1  istr ibutions  for  these  reactions  in  order 
to  describe  the  distributions  in  terms  of  a  small  set  of  para¬ 
meters,  and  to  predict  the  full  vibrotational  distribution  for 
the  reaction 

D  '  DF(v,J)  +  F  (IV) 


We  also  apply  the  model  to  the  experimental  distributions 
obtained  by  Anlauf,  et.  al,  (reference  8)  for  the  reactions 
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The  model  is  then  usee  to  predict  the  product  c  i  ?  tr  i  bu  t  lor.  for 

D  +■  Br^  -  DBr  (V,  Jl  Br  (VTII 

This  reaction  product  distribution  has  also  not  been  studied  ex¬ 
perimentally  . 

The  model  incorporates  well-docu:r.ented  reaction  complex 
dynamics  characteristic  of  Reactions  (I) -(VIII)  along  with  ideas 
similar  to  statistical  collision  theory  as  reviewed  in  reference 
9.  The  surprisal  rechnique  is  used  to  make  the  product  distri¬ 
butions  more  tractable.  We  now  review  the  surprisal  approach 
and  apply  it  to  each  reaction. 


II.  INFO  RblA  T 1 0 :  J  C  0  ’ ;  T : 


OF  AN  EXPERiy.ENTAL 


'ION 


FroTi  development  in  reference  5  .ve  consider  the 

distribution  of  outcomes  in  an  experim.ent  which  consists  of  a 
large  number  of  repetitions,  N,  of  the  same  event,  with  n  possible 
outcomes.  The  probability  of  the  ith  outcome  is  defined  as 

=  Nj^/N,  the  fraction  of  times  that  we  observe  the  ith  outcome. 
The  infcrm.ation  content  of  the  distribution  is  then  given  as: 

P  . 


The  information  attains  its  smallest  value  (zero)  when 
P^  =  1/n  so  that  all  the  possible  outcomes  arise  with  equal 
probabi 1 i ty . 

The  entropy  of  tne  distribution  is  defined  by 

H  =  -  '  P.  .'n  r  .  (2) 

",  1-  1 


where  the  thermoaynamic  entropy  S  =  RH ,  and  R  is  the  gas  constant. 
Then 

I  -  ■n(n)  -  H  (3) 


Both  the  information  and  entropy  are  ncnnivia  1 1 -.-e  numbers. 

The  prior  distribution,  P°,  of  p'rc.luct  states  (outcomes) 
is  given  as  the  one  which  corresponds  to  T  -  0.  Thus,  for  simple 
microcanon ical  systems, 

„o  1 


1 

n 


(4a) 


or , 


if  the  states  are  grouped  closely  enough  to  dcrir.e  a  density 
of  states  function  .'(E)  ,  of  the  total  energy  E,  we  have, 

F°  =  1/t  (E)  ( 4b) 

Using  (4)  in  (1),  and  further  defining  the  "surprisal" 


I  .  =  -f n{P . /P°} 
1  11 


i  5 ; 


we  have 

I  =  (6) 

We  can  see  that  the  information  is  the  negative  average  value 
of  the  surprisal,  a  guantity  which  from  information  theory 
is  the  difference  of  ihe  self-information  of  the  exper im.ei  ta  1 
distribution  (-inP.)  and  the  self  -  inf  orma  t  ion  of  ‘.he  .microcanonica^ 
distribution  {-;nP7'. 

Since  P.  reoresents  some  combination  of  statistics  and 
1 

dynamics  in  an  expcri.'ient ,  and  P^  is  a  purely  statistical  quantity, 
the  form  of  the  surprisal  should  then  give  us  relevant  dynamical 
information . 


III.  SURPRISAL  ANALYSIS 


To  characterise  tiie  pumping  distributions  we  now  calculate 

the  form  of  the  surprisals.  The  experimental  rates  are  those  of 

references  6-7  and  vie  use  the  RRHO  approximation  to  generate  simple 

baseline  statistical  lates  P°.  We  em.ploy  the  dimensionless  energies 

f  ,  f  ,  and  f  which  correspond  to  the  fractional  product  energies 
R  V  I 

in  rotation,  vibration,  and  translation,  respect ively .  Thus  we 
must  always  have 


'  V 


f„  =  1 


(7) 


The  (f  }  are  ca leu Litod  according  to  the  total 
cons  tra i nt 

E,  .  =  -  +  E  4-  f  RT 

to-'.;  a  2 


:?ncrgies  with  the 
(8) 


where  the  hr.at  of  :  ■  •act  ion,  -YH^,  is  added  to  the  relative 

reactant  trans  lat  ;  .u'.a  1  energy,  E  +  ^-RT ,  plus  an  additional  RT 

3  2. 

for  the  int(?rnal  '■.’.’■rgy  of  the  incident  diatomic. 


Rather  than  workinq  directly  '.vith  the  rotationux  surtrisals, 
we  now  compute  and  analyze  the  translational  surprisals,  therebv 
considering  the  vibrational  and  rorational  cegrcGS  of  freedom  as 
components  of  the  internal  molecular  motion  such  that 


"  V  ^  ^  -  ^inf 


(9) 


The  translational  surprisals  for  different  viorational  levels 
are  compared  by  normalizing  P  and  F°  as  though  each  product 
vibrational  level  comprised  a  separate  experiment. 

Thus,  we  constrain 

1 


P  (  f  „  V  )  = 


:p^(f.^;v)  -  1 

where  the  probabilities  are  given 
energy,  given  a  vibrational  level 
computed  according  to  reference  5 
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where  £  is  determined  to  satisfy  the  normal i za t icn  by  selecting 
those  values  of  f..,  which  correspond  to  experim.entallv  observed 

i  “  “ 

J-levels.  VJith  this  normalization  method,  we  can  analyze  the 
deviation  from  the  statistical  rate  within  each  ■.•-level  without 
concern  for  the  weighting  effects  of  the  vibrational  distribution. 

The  translational  surprisals  for  Reactions  II)  ,  (II)  ,  (III)  , 
(V),  (VI),  and  (VII)  are  shown  in  Fi'gures  1-2.  It  is  important 

to  note  that  the  curves  are  most  narrow  for  the-  highest  ob¬ 
served  vibrational  levels  and  broaden  for  lower  levels.  In 
general,  the  surprisals  achieve  maximum  width  in  the  intermediate 
vibrational  levels.  We  also  note  the  lack  of  s-.’-metry  in  the 
surprisals.  Particularly  for  Reactions  (III)  and  (VII)  the 
surprisals  skew  toward  the  high  rotational  (lo'w  translational) 
levels,  with  respect  to  the  .most  probable  J-le-.n:-!,  denoted  J. 

We  use  vJ  for  the  rotational  level  at  which  the  .surprisal,  I, 
is  a  minimum.  The  value  of  J  is  usually,  but  not  nc'cessar  i  ly , 
equal  to  that  rotational  level  within  a  v  ib:  1 1  i  or;  a  1  manifold 
with  the  greatest  population. 

We  quantify  the  full-width  of  each  suiguisai  by  measuring 


arbitrarily  at  I  =  ln2  ,  whore  P  -  P^/2.  V.'c  see  in  Table  I 
that  for  Reactions  (I)  and  (III)  the  full-widths  are  greatest 
in  the  middle  vibrational  levels.  According  to  the  model 
described  in  the  next  section,  the  Reaction  (II)  maximum 
full  width  should  also  occur  in  a  middle  vibrational  level; 
we  justify  this  identification  later. 

IV.  A  SEMI-EMPIRICAL  MODEL  TO  CORRELATE  ROTATIONAL  SURPRISAL 

FUNCTIONS . 

As  a  result  of  the  normalization,  all  surprisals  for  a 
given  reaction  have  roughly  equivalent  peak  m.aqr. i tudes .  The 
significant  variations  are  in  the  full-widths  of  the  curves 
at  some  arbitrary  value  of  I,  and  in  the  degree  of  skewness. 

We  can  directlv  relate  the  surprisal  widths  to  the  entropy  of 
the  distribution.  Clearly,  as  the  width  shrinxs  to  a  single 
J-level,  we  have  a  minimum  of  entropy  (or  a  maxim.umi  of  information) 
i.e.,  we  can  make  the  best  possible  rotational  le\’el  prediction 
for  an  experiment.  As  the  surprisal  broadens  ultimately  to  a 
horizontal  line  such  that  I  -  0,  we  have  P  =  P°  and  thus  a 
maximum  of  entropy  since  the  result  is  simply  microcanonical . 

In  a  future  paper,  we  .quantify  the  relationship  of  the  surprisal 
widths  to  the  classical  "information"  as  defined  in  Eq .  6,  and 
show  how  the  choice  of  rate  normalization  is  crucial  to  the 
use  of  information  as  a  predictive  tool. 

We  now  construct  a  model  to  interpret  the  surprisal  results 
based  on  reaction  complex  dynamics  as  described  in  references 
10-14.  Trajectory  calculations  for  three  body  exothermic  re- 

4- 

actions  of  the  typo  A  ^  BC  -  AB  +  C  from  reference  10  indicate 
that  attractive  potential  energy  surfaces  facilitate  multiple 
encounters  among  the  Liiroe  atoms  before  separation.  A  limit 
is  eventually  reached  where  the  dissociation  of  the  complex  is 
governed  by  statistical  considerations,  as  o-xplored  by  J.  C. 

Light  and  cc-workers  in  references  15-18.  In  cases  where  the 
potential  energy  surfaces  exhibit  a  mixture  of  attractive  and 
repulsive  character,  we  find  useful  both  statistical  and  dy¬ 
namical  concepts. 


The  net  effect  of  secondary  encounters  v;itr.ir.  the  reactio:. 
complex  is  a  broadening  of  the  energy  and  angular  uistributions , 
as  proposed  in  references  10-14.  The  two  primary  types  of  en¬ 
counters  are  termed  "clouting"  and  "clutching."  Clouting  in¬ 
volves  a  repulsive  encounter  between  A  and  C  '.-.•ith  conseauent 
reduction  in  the  angular  velocity  of  AB  .  Clutching  secondary 
encounters  are  also  known  as  "migratory"  encounters.  The 
probability  of  "migration"  increases  with  the  magnitude  of  the 
attractive  part  of  the  potential  energy  surface,  A  .  Therefore, 
for  exa.mple,  we  would  expect  m.icration  to  play  a  creater  role 
in  Reactions  (III)  and  (VII)  where  A  -45%  than  in  Reaction  (V) 
where  A  =25%  (references  7  and  19) . 

A  furtner  effect  cited  in  references  10-14  is  an  observed 
decrease  in  the  mean  [.roduct  vibrational  energy,  ' 
increasing  collision  complexity  resulting  from,  increasing  A  . 

The  result  is  conversion  of  vibrational  energy  to  rotation  and 
translation.  Product  rctationai  energy  is  particularly  enhanced 
for  migratory  encounters. 

Dramatic  migratory  effects  for  the  reactions  ii  ICl  h'Cl  + 

H  +  ClBr  *  (HCl  +•  Br,  HBr  +  Cl)  resulting  in  bimodal  rotational 
distributions  within  each  product  vibrational  level  are  described 
in  reference  13.  Applying  these  concepts  to  Reactions  (I)-(,VIII), 
we  anticipate,  and  in  fact,  observe  both  broaaened  and  skewed 
rotational  distributions  with  the  most  attractive  potential 
energy  surfaces  exhibiting  the  largest  asymm.etry  toward  high 
rotational  levels. 

A  model  which  combines  these  dynamical  ideas  witn  statistic- 
must  predict  the  vibrational  dependence  of  be  r,  the  surprisai 
full-widths  and  the  Jovgree  of  asymmetry  or  skewness.  A  schematic 
of  the  model  is  shown  in  Figure  3.  For  pedogouical  purposes 
Reaction  (I)  is  used  since  only  three  v- levels  are  invol\'ed. 

As  the  reaction  begins,  a  particular  prodv.ct  "lyirational 
energy  becomes  most  probable  depending  on  a  giv.in  trajectory 
through  the  potential  energy  surface.  This  r  a  t  lonal  energy 
may  be  quickly  transformed  into  rotational  and  o.sl  ational 
energy  via  clouting  or  clutching  secondary  .-ruvsun t  (--rs  followed 
by  repulsive  energy  reli^ase. 


In  order  to  modei  these  interocr icr.s  wc  ejr ^nau  th-. 
encrqy  redistribution  in  vibration,  rotation,  and  t lan.s xat ion 
(V,  R,  and  T)  can  be  described  by  a  redistr  it'ur  :  or.  o:  i’  and  F 


states  within  the  fully  formed  moleculG--HF  for  example  in 
Reactions  (I)  and  (III) .  Use  of  the  product  V  and  R  states  is 
especially  reasonable  in  Reactions  ( 1 1 1 ( VI 1 1  ;  in  which  ‘he 
"light-atom  anomaly"  results  in  an  "AB"  bond  being  effectively 
formed  before  the  "C"  a xOm  becomes  involved.  This  allows  use 
of  the  unperturbed  product  states  as  an  approx  im.a 1 1 or.  to  the 
time-dependent  reaction  complex  wave  fur.cticns. 

For  the  st  '/cifio  case  shown  in  F:  :ure  3  we  .-ce  that  the 
v  =  l  population  will  result  partly  from.  :'0 1 'xcu .  initially 
choosing  the  .•  =  1  product  level  and  pa;‘iy  from  m;i'-cules  choising 
levels  v  =  2  and  ■.'-3  me  undcruoing  the  enerciy  re  d  ;  s  t  r  ibu  1 1  on 
process.  The  red  istr  ;  out  i  on  path  is  taken  tr  :e  from,  vibremion 
resonantly  to  rotation  and  then  in  part  to  trmsiation.  This 
would  seem  to  follow  the  order  taken  in  the  cem.plex.  A  direct 
path  from  vibration  tc  a  mixture  of  rotation  and  translation 
was  also  tested  with  ign  i  f  icant  ly  inferior  results. 

The  contributions  from  the  higher  V-levels  must  be  weighted 
according  to  the  initial  transition  state  vibrational 
probability,  and  the  amount  of  rotational  energy,  *E, 

converted  into  product  translational  energy.  The  term  E 

V  “  H 

denotes  the  energy  resonance  among  the  selected  V-R  states. 

If  the  initial  populations  are  taken  as  grouped  wiL’nin  equally 
small,  energy  intervals.  A' ,  the  redistribution  process  makes 
these  states  available  Lu  lower  V-leveis.  The  natural  logarithm 
of  the  number  of  wei  dated  available  states  then  eorresponds  to 
the  entropy,  or  width,  of  the  given  distribution. 

We  take  P(E,,  „)  Luiual  to  the  final  vibrational  population, 
P(v)  .  We  find  thii.  ilthough  the  model  predicts  large 
difference  between  'he  ib.sol  ute  reaction  comp'lex  and  final 
vibrational  poru  1  ,i  1 1  ons  ,  the  relative  (normalised)  [copulations 
are  almost  unchanged.  We  see  in  fact  a  s.mall  relative 
population  shift  to  lower  v-levels,  but  tfie  iiffcrenccs  are 
actually  within  the  ■' xper  imental  uncertainties  of  the  vibrational 


distribu'ior  Cach  state  is  aise  ’.ve  i  ah  tea  by  .. 
which  gi'  .  the  probability  lor  rciaxaticn  :r:r  sel.stoa 
state  wo  the  n.ost  probable  rotaticiial  state,  a',  at.  ■.r.ei  r.- 
distance  : E  away.  The  adjustable  parameter  "  is  a  measure 
of  the  magnitude  of  the  energy  redistr  ibu  t  i.ur.  iiacess  far  a 
given  reaction.  It  is  analogous  to  the  reciprciai  of  the  con¬ 
stant  "C.,"  as  given  by  Polanyi  and  Woodall  in  reference  20  for 
HE  collisional  rotational  relaxation: 
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( see 


We  select  an  enercv  interval,  1 


;  luure  ,  acnrcxi- 


mately  equal  to  twice  the  smallest  surer isal  width  and  then  vary 
parameter  a.  (In  practice  the  specific  choice  of  interval  has 
little  effect  on  the  results  as  long  as  it  is  reasonably  small.', 
For  small  a,  we  find  that  high  J-levels  are  not  contributing 
enough  to  the  rotational  width,  while  at  large  .  the  ccntributio: 
is  too  great.  The  optimum  value  of  a  is  determined  by  locating 
the  minimum  in  the  standard  deviation  function  i  =  :(a)  where: 


(V,  0 


oreaictec 


( V  ,  . ; 


ocservec 


In  Reaction  (Illi  for  e.vample,  using  data  for  the  first  fi'oe 
vibrational  levels  results  in  a  pronounced  r.i.nimum  for  r;  ./ 
at  a.  =  14.  For  Reacoion  , ,  althcuch  thorn  are  wiaths  to 
correlate  at  only  v  =  I  and  v  =  2,  we  fine  an  optim.um  x  =  6. 

Here  we  used  the  v  =  u  surprisal  width  as  naif  the  energy  i;'. torva 

The  Reaction  (II/  rotational  distribution  for  v=2  was 
determined  by  Polar.yi  (reference  6)  with  considerably  less 
accuracy  than  for  the  v  =  3  and  v  =  4  levels.  Furthermore, 
the  V  =  1  distribution  was  esti.mated  since  he  observed  no 
detectable  emission  in  that  band.  If  we  apply  our  techniques  to 
Reaction  (II),  assumirvg  =  6,  we  can  then  generate  a  rotational 
distribution  for  the  first  vibrational  level  as  shown  in 
Fig.  5.  The  data  ■..•■xs  generated  using  parameters  from  the 
other  vibralional  lu-vels.  This  results  in  large  entropy 
rotational  distributions  in  the  middle  v- levels  for  Reaction 

(II)  as  we  iiavG  already  observed  for  Reactions  (I)  and 

(III) .  We  chow  th-  results  for  Reactions  (I)  ind  (II)  in 
Table  II. 

The  pururujter  " is  taken  physically  as  the  strength  of 
interaction  for  f  h.e  ceco.ndary  encounters  in  !  he  i/iodol.  I.arue 
values  of  .  'i'.'not.c  higiv'r  probabilities  :'or  ct'd  i  s  tr  ibu  t  ion 
over  s  i  gn  1  f  1  n  t  •  ala  'c  af  'F..  ac,-  ■  ^hat  x  is  dependent 

on  the  :  ct.ntial  energy  surfaces  .and  the  mcicroscopic  reaction 
initial  coi.  iitions.  Thun  x  may  well  be  ip,:  r./ximately  isotocac- 
ally  indet,  -ndent  ,  rarticularly  botwe.jn  Rt;actions  (III)  and  (IV) 


Wc  note  that 


and  bctwonn  (VII)  anc  (’■■III)  .  Wc  note  that  tnc  octirr.ur,  vaiacs 
of  both  A  and  (.  for  Peactions  (V)  and  (VI)  are  very  close,  in¬ 
dicating  tha*-  isotopic  independence  of  these  pdru;a,ters  vs  a 
reasonable  assu;nption.  Using  this  assu.Tiption  we  can  generate 
the  surprisal  widths  predicted  by  the  model  for  reactions  (IV) 
and  (Vill)  as  shown  in  Figure  6.  It  is  not  surprising  that  the 
widths  for  Reactions  (IV)  and  (VIII)  generally  exceed  those  of 
Rec  tions  (III)  and  (VII'  since  the  deuterium  product  energy 
levels  are  closer  together,  facilvtating  energy  level  mixi:,g. 

As  an  almost  independent  check  on  the  model,  we  now  predict 
the  degree  of  asymmetry  for  rotational  surprisals  in  reactions 
characterized  bv  attractive  surfaces.  We  excoct  rotational 


asymmetry  in  migratory  encounters  which  may  enhance  iiroduct 
rotation  as  opposed  to  repulsive,  clouting  .-nccunters  which 
tend  to  restrict  predict  ictational  energy.  If  we  now  compare 
the  probability  of  energy  transfer  to  a  low  J-levei  (at  say 
I  =  in2)  with  that  of  a  high  J-levei  (also  ^t  I  =  ;n2)  we 
obtain  a  qiantitative  measure  of  the  asymumecry,  i.e.,  a  com¬ 
parison  of  the  " ha  1  f-w laths "  of  the  distributio.'s  centered 
about  eacn  J. 


Therefore  consider 
which  transforms  within 
energy  E^,  or  to  a  high 
these  two  p^robao 1 1  i ties 
analogy  is  then 


an  intermediate  state  with  energy 
a  v-level  either  to  a  low  J-level  of 
J- level  of  energy  E.,.  The  ratio  of 
according  to  a  collisional  relaxation 
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(17 


If  E2  -  E.  IS  given  as  the  r-nergy  width  of  a  tr.Anslational  sur¬ 
prisal,  then  we  can  pr-;'dict  the  ratio  of  "  na  1  f -w  i  dths  "  as 


1  -(E^-E,)/ikT 

•  ^  n 

whore : 

hj  =  f^(I— ;n2  at  low  J)  -  f^(J) 

■'2  “  (I  -^2  at  higih  ,1) 


(IS 


(  1  9  1 ) 
(  1  ^'b) 


f  ^  ( I  -  ■  nf;  at  low  J) 
f„(I=.n2  at  hiah  J) 


(  19c; 


(  19d; 


A  4.  ‘ 

*^1  “2 


=  E2  -  El 


(19e) 


A  comparison  of  Equation  18  with  the  actual  half-widths  for 
Reaction  (III)  is  shown  in  Figure  7.  We  see  that  for  tne  same 
values  of  a  we  obtain  good  predictions  for  both  the  half-widths 
and  the  full-widths.  The  same  technique  was  applied  to  Reaction 
(VII)  with  less  quantitative  success  due  primarily  to  the  sim.- 
plicity  behind  equations  (17) -(19)  in  estim.ating  redistribution 
probabilities.  The  .-'odel  does  exhibit  qualitative  success  in 
predicting  Reaction  iVII)  asymmetry  and  ccul.i  be  m.ade  more 


soDhi s  t ica  t  e 


that  our  inte.nt  is  t< 


rovide  a  simcle 


m.odel  for  correlating  product  distributions. 

Vve  now  develop  a  sim.ple  algorithm  to  reconerate  the 
rate  coefficients  for  Rec.ction  (III)  based  on  the  model.  Using 
the  translational  surprisals  we  note  that  most  of  the  curves 
can  be  approximated  by  a  linear  function  for  f^  <  f^  (J) 
and  a  quadratic  function  for  f^  >  f^  (J) .  So  for  large  f^, 
we  find  "A"  such  that 

I  -  I  (J)  =  A[f,j,,I  =  .n2)  -  f^(J)  ]  "  (20) 

where  I  =  «ln2  and  f,„(J)  is  the  value  of  f^  at  the  maximum  value 

r  1 

(pealc)  of  [-[(f^)].  The  term.  f.p(I  =  I.n2)  is  determined  by 
Icnowing  that 

^  c^nQ'  (2i; 


^  -  ■./  .)^T 


Combining  Equations  21  and  22  gives 


Then  for  lar 


g f , 


'./  d;T 

small  f  )  we  have 
H 


.'n2  It  low  J) 


( .J )  +  ■ 


rotational 


I(J)  IS  deterrr.inec  by  observing  that  eacr.  rrociuct 
distribution  has  total  population  roughly  proportional  to  the 
width  times  the  height  of  the  distribution.  Thus,  we  set 

Y  P  ( V ) 


P(J)  = 


(26) 


where  y  is  fixed  to  make  the  largest  value  of  P(j!  for  the 
entire  vibrotational  distribution  equal  to  unity.  Now  invertir.c 
Equation  5  and  combining  with  Equation  12  we  have  for  any  value 


of  f,^: 


where 


P{f^;  =  - (v)P°(f  )e 

1  1 


<  ( v ) 


( J  i  V ) 


with  P'  as  the  experimental  distribution  obsirv'd  in  reference  ~. 
The  last  two  formulas  are  simply  the  inversion  ef  the  norm.ai  i  cation 
method  employed  previously.  For  I(J)  we  now  have: 


I  (J)  = 


,P  (J) 


P^(J) 


=  -  i  n  ( 


•,P(v)  , 


">  Q  ^ 


For  f  <  frpid)  we  find  the  equation  of  a  line  through  the 

points  I(J}]  and  [f^(J)-h2/  I  =  (n2].  Using  the 

linear  and  quadratic  form.s  for  the  surprisa.s  as  functions  of 

f^,  we  can  now  generate  the  actual  distribution  P  -  P(f,p)  by 

using  Equation  27  for  eacn  selected  value  of  f^. 

Figure  8  compares  with  our  serr. i -emu ir ica  1  model  prediction 

for  the  Re.-iction  (III'i  race  coef f icien.ts  with  t;ie  t-xper im.enta  1 

data  in  terms  of  f  .  The  model  also  uives  oualitativo  aareemont 

K 

with  the  Reaction  ('.’IT)  -jxper  i  mental  data  il  though  the  sim.plici- 
ty  of  the  algorichm  results  in  some  cuant i tat ive  differences. 


VI.  RATE  trCEFFICIi::.  r  prediction  for  PEACTIONc 
(IV)  and  (VI  ID 

We  now  ipply  the  model  with  some  additional  assumptions  to 
Reactions  (IV)  and  (VIII)  for  which,  to  the  best  of  our  knowle-dge, 
detailed  experimental  rate  coefficients  arc  not  yet  available-. 

To  develop  the  full  v’ibrotational  d  i  st  r  i  nu  t on  v;e  must  predict 
a  vibrational  pumping  distribution,  a  mo;^t  pio'Eable  J  within 


each  v-level,  and  the  form  of  the  rotational  surprisals  for 
each  v-level.  Berry  has  reported  in  reference  21  an  isotopic 
independence  for  the  vibrational  surprisals  of  Reactions  (I) 
and  (II).  The  linearity  of  these  surprisals  indicates  that 
the  distribution  is  characterized  by  just  one  moment,  namely 
the  mean  product  vibrational  energy.  The  vibrational  surprisals 
for  Reactions  (III)  and  (VII)  are  computed  using  the  RRHO  result 

P°(f^,)  =  3  (1-f^)  (30) 

such  that 

y  =  1  (31) 

The  resultant  surprisals  are  not  entirely  linear,  but  are 
regular  enough  so  that  if  we  assume  an  isotopic  independence 
we  can  predict  I(f^)  fot  the  product  DF  and  DBr  molecular 
vibrational  levels.  Then  using  the  predicted  values  I (v)  we 
find  the  vibrational  distribution 


where  is  such  that  the  largest  value  of  P(v)  is  unity. 

This  normalization  giving  P(v)  =  1  is  the  same  used  by 
Polanyi  in  reference  7. 

There  is  also  reported  (see  reference  22)  an  isotopic  in- 
deoendence  for  f  as  a  function  of  f  ^  between  Reactions  (I)  and 
(II).  Assuming  the  same  relationship  between  Reactions  (III) 
and  (IV)  and  between  (VII)  and  (VIII)  we  can  predict  J  (v)  for 
Reactions  (IV)  and  (VTII). 

To  predict  the  rotational  surprisals  according  to  the 
model,  we  use  the  widths  predicted  in  Figure  6,  and  the 
same  prescription  as  given  in  Section  V  for  Reaction  (III)  . 

We  further  identify  that  only  one  rotational  level  beyond  the 
full-width  is  selected.  The  full  vibrotational  distributions 
for  Reactions  (IV)  and  (VIII)  thereby  generated  are  shown  in 
Figures  9  and  10. 


VII.  SUMMARY 

A  model  was  developed  which  combines  statistical  ideas 
with  well-established  reaction  dynamics.  The  model  describes 
the  features  of  experimental  rotational  distributions  for  six 
reactions  and  predicts  the  full  vibrotational  distributions  for 
two  reactions  not  yet  determined  experimentally.  Future  work 
involves  defining  the  relationships  among  arbitrary  rotational 
distributions  and  the  classical  "information"  of  a  distribution 
for  a  given  reaction.  V.'e  also  hope  to  extend  the  ideas  presente 
here  to  ccllisional  relaxation  processes. 


TABLE  I  . 


Translational  Surprisal  Full-Widths, 


Reaction  (I) 

V_ 

1 
2 
3 


0.198 

0.212 

0.045 


Reaction  (II) 


V 


1 


0.181 


0.147 


3  0.122 

4  0.054 


Reaction  (III) 


V 


1 

2 

3 

4 

5 

6 

7 

8 


0.036 
0.048 
0.065 
0.068 
0.064 
0.068 
0.  030 
0.017 


^) ,  at  I=.n  2 


TABLE  II.  Translational  Surprisal  Ful l-V.’idths , 


at  1=  ; 


Reaction  (I) 


V 

1  observed 

1  predicted  v 

) 

1 

0.198 

0.208 

2 

0 . 212 

0.212 

Reaction  (II] 


V 

1  observed 

.L  predictea 

.  .  =  6) 

0.181 

0.  117 

2 

0.147 

0.147 

3 

0.  122 

0.131 

The  model  predictions  are  compared  with  the  observed  data  fo 
Reactions  (I)  and  (II).  Note  that  the  precise  match  for  the 
largest  full-widths  is  simply  due  to  the  normalization  of 
equation  (15) . 
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Figure  1.  Surprisals  as  functions  of  translational  energy  for 
each  of  the  experimentally  observed  product  vibrational  levels 
ror  Reactions  (I)-(III).  The  surprisals  were  computed  using 
the  normalization  described  in  the  text. 

Figure  2.  Translational  surprisals  for  Reactions  (V)-(VII). 

Figure  3.  Schematic  for  the  model  predicting  surprisal  widths. 

An  energy  interval,  L' ,  is  chosen  corresponding  to  the  v  =  3 
surprisal  width  at  the  value  I  =  ln2.  All  of  the  energy  levels 

A. 

within  f  1  A,  where  f  is  the  enercv  of  the  v  -  3  pooulation 

maxim, um,  define  a  group  of  states  of  approximately  equal  energy 

E,,  This  cTrouD  of  states  includes  these  fallina  within  the 

interval  at  lower  vibrational  levels.  According  to  the  energy 

redistribution  mechanism  described  in  the  text,  each  rotational 

distribution  will  include  broadening  effects  due  to  the  selected 

states  accessible  to  a  given  vibration  level.  Thus  for  the 

reaction,  the  v  =  2  level  will  have  contributions  from  the  v  =  3 

level,  and  the  v  =  1  level  will  have  contributions  from  both 

the  V  =  3  and  v  =  2  levels.  Each  selected  rotational  level  is 

weighted  according  to  level  dege.neracy  (2J+1),  population 

and  the  rotational  energy  distance,  cE ,  from  the  level 

to  the  rotational  distribution  center  f  . 

R 

Figure  4.  Predicted  (1)  and  observed  (0)  values  of  the  trans¬ 
lational  surprisal  widths  as  functions  of  the  normalized  vi¬ 
brational  energy,  f^.  The  constant  C  is  determined  by  equating 
the  maximum  predicted  and  observed  widths  for  each  reaction. 

We  expect  and  observe  the  best  results  for  Reactions  (III)  and 
(VII)  where  many  vibrational  levels  are  involved  in  the  energy 
redistribution  mechanism  which  is  the  crux  of  the  model. 

Figure  5.  The  predicted  v  =  1  product  rotational  distribution 
for  Reaction  (II),  using  the  techniques  of  t!ie  model. 

Figure  6.  The  predicted  translational  surprisal  full-v.’idths 
for  Reactions  (IV;  and  (VIII) .  The  data  is  generated  assuming 
isotopic  independence  in  the  para.mot  irs  i,  C,  and  f^^  (  -  between 


Reactions  (III)  and  (IV)  and  between  (VII)  and  (VIII). 

Figure  7.  The  ratio  of  predicted  surprisal  "half-widths"  (1) 
compared  with  observed  data  (0)  for  Reaction  (III).  Data  for 
Reaction  (VII)  shows  qualitative  agreement,  although  quantita¬ 
tive  discrepancies  arise  due  to  the  simplicity  of  characterizing 
the  rotational  energy  distance,  6E,  as  the  difference  of  each 
selected  level  energy  and  f  .  Constructing  a  model  using  a 
continuum  value  for  6E  as  the  difference  between  each  selected 
level  energy  and  weighted  values  of  the  product  rotational  dis¬ 
tribution  energies  might  improve  quantitative  results,  but  would 
also  involve  unwarranted  com.plexity  in  an  otherwise  simple  model. 

Figure  8,  Predicted  (1)  and  experimental  (0)  product  rotational 
distributions  for  Reaction  (III).  Note  that  the  results  are 
best  in  v-levels  where  the  model  accurately  predicts  the  sur¬ 
prisal  full-widths.  Reaction  (VII)  data  exhibits  qualitative 
agreement  although  quantitative  discrepancies  arise  due  to  the 
simplicity  of  the  model  techniques,  as  alluded  to  in  the  caption 
to  Figure  7. 

Figure  9.  Predicted  product  distributions  for  Reaction  (IV). 
Figure  10.  Predicted  product  distributions  for  Reaction  (VIII). 
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